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ABSTRACT 

We calculate the effect of illumination of dwarf nova accretion discs by radiation from 
a hot, central, white dwarf. We show that only for very hot white dwarfs (Teff w 
40000 K) the inner region of quiescent dwarf nova discs are partially depleted so 
that the delay between the rise to outb urst o f the optical and UV fluxes would be 
increased as suggested recently by King ( 1997 ). This depletion, however, must create 
several small outbursts between main outbursts, contrary to observations. Lower white 
dwarf temperatures may cause the outburts to be of the ‘inside-out’ type removing 
the UV delay. We conclude that white dwarf irradiation of dwarf nova discs is not 
very efficient for example because the UV radiation from the hot white dwarf does not 
penetrate deep enough in the disc atmosphere. The total ablation of the inner disc 
by e.g. evaporation (possibly related to illumination) appears to be a very promising 
possibility, accounting for both the EUV delay and the general lightcurves properties. 

Key words: accretion, accretion discs - instabilities - novae, cataclysmic variables 
- binaries : close 


1 INTRODUCTION 

Dwarf novae (DN) are cataclysmic variables which exhibit 
large amplitude outbursts that are commonly attributed to 
a thermal/viscous instability occuring in the accretion disc 
around the white dwarf (pr imary ) component of the binary 
system (see e.g. Cannizzo (1993) for a review). In this pic¬ 
ture, a stable disc has a bimodal behaviour: it can be either 
in a hot, ionized, state, corresponding to large mass transfer 
rates, or in a cool, neutral state, with a very low mass trans¬ 
fer rate. Systems in which the mass transfer rate from the 


secondf ry is in the intermediate range cannot be steady. The 
disc becomes unstable when partial ionization of hydrogen 
sets in; this may occur either in the innermost parts of the 
disc, in which case the outburst is of the inside-out type, or 
in the outer disc (outside-in outburst). In both cases, a heat 
front forms and rapidly propagates across the disc, bring¬ 
ing it entirely (or almost entirely) to a hot state, in which 
the local mass transfer rate is larger that the rate at which 
matter is provided by the secondary. The disc thus empties 
until a cooling wave starts from the outer edge of the disc, 
and brings it entirely to the cool (and unsteady) state. 

Whereas the disc instability picture accounts for the 
general observational properties of DNs, it still has to face a 
number of contradictions, both internal and external. Mod- 
els suffer from a number of drawbacks and incon sistencies 
(Lasota fc Hameury, 199^; Hameury et ah, 1998). The de¬ 


scription of viscosity which drives accretion onto the white 
dwarf is purely phenomenological (and often ad hoc), so that 
it may be difficult to determine if failures to reproduce some 
observational properties of dwarf nova outburst are due just 
to the approximate treatment of the disc structure, or to 
some more fundamental physical problems (see e.g. Gam- 
mie & Menou 1998). 

During some DN outbursts one observes a so called ‘UV- 
delay’ - the rise to outburst in the optical wavel ength pre¬ 
cedes the UV rise by as much as 0.5 to 0.75 day (Hassall et 


al., 1983; Warner, 1995); the d elay between th e optical and 
EUV can be as long as 1 day (Mauche, 1996). Such delays 


are observed during outburts which start at the outer disc 
regions - the ‘outside-in’ outbursts and are usually assumed 
to be the propagation time of a heating front across the disc. 
Whereas this interpretation can be questioned in the case of 
1000 - 2000 A UV, since t he spectral mode lling of the disc 
must be accurately done (Cannizzo, 1998), the EUV rise 


clearly marks the arrival of the heating wave at the white 
dwarf surface. The presence of the UV-delay has been a se¬ 
rious difficulty for all disc models, since according to them, 
the heating front propagates too rapidly across the disc - the 
disc becomes hot too fast for an appreciable delay between 
UV and optical to appear (see Fig. ^). Attempts to solve 
this problem by modifying the local physical proper ties o f 
the disc (e.g. increasing the thermal scale - Mineshige (198E) 
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- or involving convection -Duschl (1989)) have not been suc¬ 


cessful. A mu ch more promising possibility was suggested by 
Meyer (199C): he proposed that the innermost regions of the 
disc are ablated by a coronal siphon flow in quiescence, so 
that UV emission, originating from the immediate vicinity 
of the white dwarf, can occur only after the central disc has 
been reconstructed, which is happening in a viscous time 
(see Hameury et ah, 1997b for the description of the same 
phenomenon in low mass X-ray binary systems). 

Although the inability of DN models to reproduce ob¬ 
served UV-delays has been attributed to the simplified way 
disc radiation is treated in time-dependent calculations (see 
e.g. Cannizzo 1998) the EUV delay represents the real 
heating-front propagation time because EUV radiation is 
not emitted by the disc but by a boundary layer between 
the disc and the white dwarf. 

The precise physical mechanism causing the truncation 
of the disc is still a matter of debate. In addition to the 
evaporation of the inner dis c by siphon flow describ ed by 
Meyer & Meyer-Hofmeister (1994), Livio & Pringle (1992) 
suggested that the inner hole is due to the presence of a weak 
(^ 10^ G) magnetic field. It is interesting that Patterson et 
al., (1998) found that the white dwarf in WZ Sge may have 
a magnetic field (1 — 5) x 10^ G and that the presence of 
a truncated disc in this system had been proposed before 
by Lasota, Hameury & Hure (1995) and Hameury, Lasota 
& Hure (1997a). 


In a recent paper. King (1997) argues that irradiation 
by the white dwarf also l eads to the truncation of the disc. 
According to King (1997) disc regions close to a hot enough 


white dwarf cannot be brought to a cool state at the end of a 
DN outburst, so that matter is transferred rapidly onto the 
white dwarf, leading to the formation of a very low density, 
hot flow. 

The effect of disc irradiation in the context of soft X-ray 
t ransi ents was studied by Tuchman, Mineshige & Wheeler 
(1990), and may pl ay an important ro l e in low-mass X-ra y 
binaries in general (van Paradijs, 1996; Dubus et ah, 1998). 


In this paper, we study the effect of the illumination of the 
disc b y the white dwarf. We show that, as proposed by King 
(1997), the inermost disc regions are depleted, but contrary 
to his hypothesis they still remain optically thick. The main 
effect of this illumination is to render the transition region 
between the hot and cold regions extremely unstable. As a 
consequence, many small outbursts are triggered, outside-in 
outbursts are difficult to obtain, and most often, the para¬ 
doxical result of illumi nation is a suppression of the UV 
delay. In Dubus et al. (1998) we present a detailed analy¬ 
sis of the case in which the illuminating flux originates from 
the disc itself, and is therefore coupled to the mass accretion 
rate onto the compact object; in the present study we find 
that adding irradiation due to accretion does not modify our 
main concfusions. 


2 MODIFICATION OF THE S - Tefr CURVES 
BY IRRADIATION 

The main effect of disc irradiation from the outside is to 
change the surface temperature. A fraction 1 — /3 of the inci¬ 
dent flux is absorbed in optically thick regions, thermalized 
and reemitted as photospheric radiation. The remaining is 



E (g cm ) 

Figure 1. S — Tefr curves for a = a{T,r) (see text), a 0.6 Mq, 
8.5 X 10* cm radius white dwarf, with (1 — = 25,000 K, 

for radii 0.85, 1.09, 1.29, 1.39, 1.54, 1.71, 1.90, 2.11, 2.42 xlO® 
cm from left to right. Note that Tefr accounts only for viscous 
dissipation and differs from the actual surface temperature. 


either scattered, or absorbed in optically thin regio ns, possi¬ 
bly forming a warm corona (tdan & Shaviv, 1996). Because 


the white dwarf temperat ure is not very different from that 
of the disc surface. King (1997) assumed that the albedo /3 
is not close to unity - in this case effects of irradiation on 
the inner disc structure can indeed be important. 

In the case of a geometrically thin disc surrounding a 
white dwarf whose radius is larger than the disc thickness, 
the illumination temperature Tm defined as 


(UTlti = Fill 


( 1 ) 


where Fn is the radiative flux illuminating each side of the 
disc is given by (Friedjung 1985; Smak 1989; Hubeny 1991): 

Till = (1 - /3)r*-[arcsinp - p(l - (2) 

TT 

where p — Rt/r, R, and T* are the white dwarf radius and 
temperature, and r is the radial coordinate in the disc. 

We have solved the disc vertical structure in the opti- 
cally thick approximation, as described in Hameury et al. 
( 1998| ), having modified the disc outer boundary condition 
according to: 

= Fvisc + aFti (3) 

where Ts is the disc photospheric temperature and Fvisc is 
the energy flux due to viscous dissipation. Examples of the 
resulting E — Ten curves are given in Fig. 0 where = 
Fvisc (see e.g. Dubus et ah, (1998)), which is proportional 
to the mass transfer rate in steady state. We have assumed 
a bimodal behaviour of the Shakura-Sunyaev parameter a, 
taken as: 


log(a) = log(acoid) 


[log(ahot) - log(acoid)l/ 
/ 2.5 X 10“* K ^ ® 

(F? + 5Ti2 )i/2 


(4) 


where Tc is the central disc temperature, and acoid and Ohot 
are constants, here taken to be 0.04 and 0.2 respectively. 
Note the presence of a Tm term, that did not appear in 
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Figure 2. Maximum value of E on the lower stable branch of 
the S — Teff curve. The primary is a 0.6 Mq, 8.5 X 10® cm white 
dwarf, with effective temperature 25,000 K (lower solid curve), 
or 15,000 K (upper solid curve). The unilluminated case (dashed 
curve) is also shown for comparison. 
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Figure 3. Visual light curves of the accretion disc with a 0.6 
Mq, 8.5 X 10® cm radius white dwarf, M = 10^® g s“l with 
(1 — = 25,000 K (upper panel), 15,000 K (intermediate 

panel), and no illumination (lower panel). The contribution of the 
primary or the secondary is not taken into account 


Hameury et al. ( 1998| ), and which is needed to have a = ahot 
on the hot branch for moderate Tc. Fig. shows the the 
‘effective S-curves’ used in time-dependent calculations and 
not the constant-a S-curves that a re us ually found in the 
literature (see also Hameury et al. (1998)) 

As can be seen, regions in the immediate vicinity of 
the white dwarf are stabilized by irradiation from the white 
dwarf when the central disc temperature is pushed above the 
hydrogen ionization range. The condit ion for this to happen 
is close to that given by King (1997), i.e. Ts > 6500 K as 
in unilluminated discs, but with Tg now given by Eq. ^ this 
condition slightly overestimates the effect of illumination, 
as the relation between Tg and the central disc temperature 
is altered. This is however not the sole effect of irradiation. 
Figure ^ shows the corresponding maximum value of surface 
density on the cool branch, Emax, as a function of radius, 
Q being given by (^). On Fig. 2 the curves describing maxi¬ 
mum allowed S for cold, stable, disc equilibria end at radii 
where the disc become stable (there are no more ‘S-curves’ 
for smaller radii). At larger radii, however, irradiation has a 
strong destabilising effect as seen in Fig. 2 where Emax is re¬ 
duced by more than an order of magnitude compared to the 
non-irradiated case. In these regions, irradiation increases 
the midplane temperature, bringing it close to the hydrogen 
partial ionization regime. Equivalently, the maximum mass 
transfer rate on the lower branch is severely reduced. This, 
as we will see in the next section, has an important influ¬ 
ence on the properties of dwarf nova outbursts when they 
are affected by irradiation from the white dwarf. 


3 LIGHT CURVES AND DISC IRRADIATION 


We calculated a grid of vertical structures modified by irra¬ 
diation. Such a grid provides the cooling term as a function 
of the central temperature and surface density and we we 
used i t in the numerical code described in Hameury et al. 
( 1998 ) to determine the time evolution of the outbursts. Fig¬ 


ures H and m show the general light curve of a system with 
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Figure 4. Same as Fig. ^ but for M = 10^^ g s Note that 
the scale of the time axis are quite different 


a 0.6 Mq, 8.5 x 10® cm white dwarf surrounded by a disc 
with an average outer radius 2 x 10^® cm, fed at a rate of 
10^® and 10*^^ g s“^ respectively, for three different values of 
(1 — 0, 15,000 and 25,000 K. In the unilluminated 

case, inside-out outbursts are obtained at low mass transfer 
rates, whereas one gets outside-in outbursts for M = 10^^ g 
s“^. 

The first visible effect of illumination is the appearance 
of intermediate, smaller outbursts that are always of the 
inside-out type; these appear because of the destabilization 
of the inner edge of the disc by irradiation as seen on Fig. 
^ They are at least one magnitude fainter than the ma¬ 
jor bursts, but are still detectable. During these small out¬ 
bursts, the heat front is not able to propagate across the 
entire disc, but dies out at rather small radii (typically less 
that one half of the disc radius). The characteristics of the 
main outbursts are not drastically changed, with one excep¬ 
tion. With M = 10^^ g s“^ and no illumination, the out- 
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r (lO^^'cm) 

Figure 5. Radial profiles of temperature T (in 10® K), surface 
density E and critical surface density Smax (in g cm“^ and local 
mass transfer rate (in units of 10^® g s“^ just prior to one of the 
major outbursts shown in Fig. |^for (1 — = 25,000 K. 



time (hours) 

Figure 6. Time dependence of the V magnitude (left two curves) 
and the ’’accretion rate magnitude” (right two curves, see text) 
at the onset of an outside-in outburst under the effect of illu¬ 
mination (solid curve, (1 — /3)^/‘^r* = 25, 000 K) and without 
illumination (dashed curve). The primary mass is 0.6 Mq and 
the mass transfer rate is 10^^ g s“^. 


bursts are of the outside-in type, as noted above. With T* 
= 25,000 K, the major outbursts are still of the outside-in 
type, while the smaller ones are inside-out. There is a transi¬ 
tion regime when T* = 15,000 K where only large inside-out 
outbursts are seen. For this intermediate regime, the surface 
density is everywhere close to Smax in the cold state (as in 
the non-illuminated case) but illumination is strong enough 
to destabilize regions close to the white dwarf. Heat fronts 
created in the inner regions of the disc can thus propagate 
throughout the disc and not only in a limited region. In 
other words, for these parameters the ‘small’ inside-out out¬ 
bursts become major outbursts. These light curves are very 
similar to those obtained with higher primary masses, and 
hence much shorter inner disc radii, in which small inside- 


out outbursts are also seen (Cannizzo, 1993; Hameury et al. 

199C] - - 

A second effect is the reduction of the outbursts am¬ 
plitude. They are only slightly fainter at maximum, but the 
main contribution to the reduced amplitude comes from the 
increased brightness during quiescence, that results from the 
brightening of the inner parts of the disc due to the repro¬ 
cessing of light from the wh ite dw arf. 

As suggested by King (1997), illumination does reduce 
the surface density in the inner parts of the disc. Fig. ^ 
shows the radial disc structure in a system with a 0.6 M©, 
8.5 X 10® cm white dwarf, M = 10^® g s“^, and (1 — 

= 25,000 K just prior to a major outburst. E is much re¬ 
duced in the region where illumination is important, i.e. up 
to about 2 white dwarf radii, but not to the point mak¬ 
ing the disc optically thin. For the value of M considered 
here, the outburst will be of the inside-out type, as can be 
seen from the relative position of the E and Emax curves. For 
higher mass transfer rates, outside-in outbursts are possible. 
Fig. 1^ shows the time behaviour of the onset of the bursts 
of both the disc visual magnitude, and the “accretion rate 
magnitude”, defined as = 27 — logMwd, where M„d is 
the accretion rate onto the white dwarf in g s“^, which rep¬ 
resents the EUV flux. As expected, the delay between the 
start of the rise of the disc brightness and the rise of mass 


accretion rate onto the white dwarf increases when illumi¬ 
nation is important. The increase is however not very large; 
for the parameters assumed here, the total delay amounts 
to about 5.3 hr, of which about 1.5 hr is attributable to 
the depletion of the inner disc; longer values are possible for 
higher white dwarf temperatures, or lower values of a. 

It is also interesting to note that illumination makes the 
rise of the mass accretion rate onto the white dwarf much 
less sudden than in the unilluminated case. This is due to the 
fact that close to the white dwarf, the heat front propagates 
now in a hot disc, with a constant value of a. This reduces 
the non-linearities which are responsible for the presence of 
sharply defined heat fronts. 


4 THE CASE OF SS CYG 


SS Cyg is a well st udied system in which a 1 day EUV delay 
has been observed (Mauche, 1996) that contains a hot, mas¬ 
sive white dwarf, whose temperature i s estimated to be in 
the range 34 - 40,000 K (Warner, 1995 ) for a mass of about 
1.2 Mq (Ritter & Kolb, 199S). In such a system, the effi¬ 


ciency of accretion is large, and illumination by the disc itself 
and by the boundary layer between the disc and the surface 
of the white dwarf could be important. The geometrical ex¬ 
tent of the hottest parts of the system is not well known; as 
a simplifying assumption, and in order to maximize the ef¬ 
fect of accretion related irradiation, one can assnme that the 
source of illumination is the white dwarf surface, so that Eq. 
(^) still applies if the constant temperature T* is replaced 
by Ts given by: 


^4 ^4 , GMiM 

aT, gT, + ^^^3 


(5) 


where Mi is the primary mass. The cooling term T^s has 
been calculated for several illumination temperatures, and 
a cubic splin e inte rpolation is performed, as described in 
Dubus et al. ( 1998 ). 
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Figure 7. Calculated visual magnitude light curves for parame¬ 
ters appropriate to SS Cyg. In the upper panel, the white dwarf 
temperature is kept constant, at a value of 30,000 K; in the inter¬ 
mediate panel, illumination by the boundary layer is also included 
(see text), and in the lower panel, we that the inner parts of the 
disc are evaporated, and we do not consider the effect of illumi¬ 
nation on the outer disc. 



Figure 8. Visual (left three curves, almost superimposed) and 
“accretion rate” (right three curves) magnitudes in the case of SS 
Cyg. The solid curve is obtained when assuming evaporation of 
the inner disc and no illumination; the dashed line corresponds 
to illumination by both the white dwarf and the boundary layer, 
and the dotted line to illumination by the white dwarf only (see 
text for parameters) 


We show in Fig. Q the calculated light curves for pa¬ 
rameters appropriate to the case of SS Cyg. We have taken 
ac = 0.02, ah = 0.10, M = 1.510^'^ g s"\ T. = 30,000 K, 
/3 = 0.5, Ml = 1.2 Mq and i?* = 5 x 10® cm. We have also 
calculated light curves for T* = 40,000 K, which are quite 
similar to those shown here. A zoom of the outburst rise 
is shown in Fig. ^ For comparison, we have also shown a 
case in which illumination by the boundary layer has been 
neglected, and a case in which no illumination is assumed, 
but matter from the disc is evaporated at a rate: 


by the boundary layer also produce light curves which show 
many small outbursts in between major ones, and cannot 
therefore account for the observed light curves. The situation 
is slightly better when the boundary layer is included, but 
is far from being satisfactory. Similarly, the predicted EUV 
delays are too short, since too small a fraction of the disc is 
affected by illumination. 

By contrast, both the light curves and the EUV delay 
obtained when one assumes that the disc is truncated are 
in good agreement with observations. The very significant 
difference between the case of a hole in the disc and a de¬ 
pletion of the inner regions is precisely the existence, in the 
latter case, of a hot, optically thick, central region that is 
able to interact with the outer disc and destabilize it. 

One should finally note that the alternation of various 
types of outbursts in SS Cyg cannot be explained if no pa¬ 
rameter (e.g. M) is allowed to vary with time. Cannizzo 
( 1993|) found that regular sequences of short and long out- 


Eev — 


3 X 10” 


(l-b250i?2J2 


-1 -2 
g s cm 


( 6 ) 


bursts were a natural outcome of the disc instability model 
for low values of the disc inner radius; one should however 
keep in mind that all his outbursts were of the inside-out 
type (because the outer disc radius was fixed at a given 
distance), which is very difficult to reconcile with the ex¬ 
istence of UV and EUV delays. We do also obtain such 
sequences, even with alternating inside-out and outside-in 
outbursts (see e.g. fig. , but the small outbursts are much 
fainter than the major ones, in contradiction with obser¬ 
vations. The long term light curve of SS Cyg shows large 
variations of the m ass transfer rate from the secondary on 
timescales of years ( [Cannizzo fc Mattel, 1992 ); these varia¬ 
tions must be taken into account when trying to interpret 
the alternation of outside-in and inside-out outbursts, and 
we know that changes in M induce changes of the outburst 
type (see Eigs. H andW). 


5 CONCLUSION 

We have shown that efficient illumination of the accretion 
disc by a hot white dwarf affects the dwarf novae outbursts 
in several ways: 

(i) it causes the appearance of several small outbursts in 
between the major ones. This conclusion is firm and does 
not depend on the assumed viscosity, but merely reflects the 
unavoidable presence in quiescence of an unstable transition 
region between the hot and cold regimes. 
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(ii) for moderate white dwarf temperatures and mass 
transfer rates, outside-in outbursts are replaced by inside- 
out outbursts 

(iii) when outside-in outbursts are present, the UV delay 
is increased. This happens only for T*(l — much above 
15,000 K. 


The proposition of King (1997) that illumination is respon¬ 
sible for the observable UV delay is therefore difficult to 
maintain since the UV delay is increased only for very hot 
white dwarfs, but then it implies the presence of small out¬ 
bursts in between major ones that are not observed. 

However, as pointed out by King (1997), such hot white 
dwarfs are present in some DNs, and their radiation does il¬ 
luminate the innermost parts of the accretion disc. We sug¬ 
gest that the fraction 1 — of the illumination flux that 
affects the vertical structure of the disc might be smaller 
than previously estimated, as (i) the opacities are a very 
sensitive function of temperature when hydrogen is partly 
ionized (factors of the order of 2 in the Rosseland opac¬ 
ity for a temperature difference of only 20%), and (ii) the 
grazing angle is quite small; both effects reduce penetra¬ 
tion depth of the UV photons, which could be absorbed in 
regions which are not optically thick and would not mod¬ 
ify the photospheric boundary condition but instead lead to 
the formation of a corona above the cool disc. The determi¬ 
nation of the albedo is, however, much more complex than 
in the case of illumination by hard X-rays, where electron 
scattering contributes significantly to the opacities. 

A much more promising explanation is that the inner 
parts of the disc that would be affected by illumination do 
not exist, precisely because of the formation of a corona 
that evaporates the central regions of the disc. Irradiation 
would indeed accelerate the evaporation process itself, but 
the coupling between evaporation and illumination remains 
to be understood. 

One should finally keep in mind the possibility that the 
disc instability model does not apply to dwa rf no va out¬ 
bursts, but then the scenario proposed by King (1997) would 
not apply anyway. 

When this work was almos t com pleted we learned of 
similar study by Stehle & King (1998). 
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